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ABSTRACT 

The elastic buckling load of simply supported rectangular orthotropic plates subjected to a 
second degree parabolic variation of axial stresses in the longitudinal direction is calculated 
using analytical methods. The variation of axial stresses is equilibrated by nonuniform shear 
stresses along the plate edges and transverse normal stresses. The influence of the aspect ratio 
is examined, and the results are compared with plates subjected to uniform axial stresses. 


INTRODUCTION 

In most studies of stability of plates, the axial stress has been taken as uniform compression 
throughout flat rectangular plates (Timoshenko and Gere, 1961; Allen and Bulson, 1980; 
Whitney, 1987). Buckling of isotropic plates under a compressive stress that varies linearly 
from one loaded edge to the other has been studied by Libove, Ferdman and Reusch (1949). 
Cases of practical interest exist, however, in which the axial stress is not uniform but varies 
from tension at both loaded edges to compression in the middle. An example is the stability of 
the crown of the hat stiffened panel, a candidate configuration of the upper and lower skin of 
the Blended Wing Body (BWB) Aircraft The BWB Aircraft is an advanced long-range ultra- 
high-capacity airliner with the principal feature being the wide double-deck body which is 
blended into the wing (Popular Science Magazine, 1995). 

In the present paper, analytical methods are used to investigate the local stability of the hat 
crown plate in order to minimiz e its weight while op timizing its buckling strength. A varying 
tension-compression-tension stress is induced in the crown of the stiffeners. The axial stresses 
vary longitudinally due to bending, and the change in stresses is equilibrated by nonuniform 
shear stresses along the plate edges and transverse normal stresses as shown in Fig. 1. 


ANALYTICAL MODEL 

The distribution of the axial stresses is given by 
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and the expressions for the shear and transverse normal stresses T^Oc.y) and <7 y (y), are 
obtained by integrating the plane stress equilibrium equations 


and are found to be 
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Figure 1. Plate subjected to tensile/co mpression/tensile longitudinal stresses o^ix ) , 
equilibrated by shear stresses r^loc.y) and transverse stresses Oyy(y) 
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The differential equation of the deflection surface of the orthotropic plate can be written as 
follows: 
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where D tJ , w, N x , N^, and N y are the bending stiffnesses of the orthotropic plate (Jones, 

1975), vertical displacement, axial, shear and transverse stress resultants, respectively A 
solution to Eq. (6) that satisfies the boundary conditions cannot be obtained in closed form 
Thus, approximate solutions are sought. Two different displacement fields are considered- ( 1 ) 
Double Fourier senes, using Rayleigh-Ritz method, and (2) a polynomial in the longitudinal 
direction, using r.alerkin method. A summary of the analytical models is described in the 
next sections, the influence of the aspect ratio is examined, and the results are compared with 
plates subjected to uniform axial loads. 


DOUBLE FOURIER SERIES (RAYLEIGH-RITZ METHOD) 

The total potential energy functional of the orthotropic plate is given by 




where t,a, and b represent the thickness, length and width of the plate, respectively. 

Of all conceivable buckling patterns satisfying boundary conditions, the actual buckling 
pattern is that for which the potential energy U , as given by Eq. (7), is a minimum. The 
deflection surface of the simply supported skin in bending may be represented by a Double 
Fourier series as 
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which satisfies the simply supported boundary conditions. Rayleigh-Ritz method requires that 
the coefficients be chosen such that U is a minimum. That is, they must satisfy the 
equations 
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The deflection expression (8) can be splitted into two independent subsets, one 
corresponding to even values of n and the other, to odd values of rt. The first subset 
corresponds to buckling that is antisymmetric about the line y = bf 2 and the second, to 
buckling that is symmetric about that line. The second subset gives lower buckling loads, and 
attention is hereinafter confined to it. Similar argument can be made to the buckling pattern in 
the longitudinal direction (subscript m). Therefore, the buckling pattern is chosen as 


^ ££ . mitx . nlty 
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Substituting Eqs. (1), (4-5), and (10) into the energy expression (7), while using Eq. (9), 
we obtain the following relations 
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P and represent the aspect ratio and buckling load, respectively. A nondimensional 
buckling coefficient k may be expressed as 
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Relations (11) represent a set of (A/ + \)(N + 1)/4 linear equations, the smallest eigenvalue 
of Eqs. (11) gives the buckling load. 


POLYNOMIAL (GALERKIN METHOD) 

We seek a nontrivial solution to the deflection w which satisfies equation (6) and the 
boundary conditions. For simply supported plates, the solution of Eq. (6) may be represented 

by _ 


* = I Z4nm(t m+3 +AnlZ m+2 +\t2t m+1 +An3Z m +K4Z m ‘ )sin(/l7T7j) 
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in which g - x j a \ and the coefficients ~ A m4 are determined by satisfying the 

four simply supported boundary conditions at the two loaded sides. 

The expression for the Galerkin equations is given by 
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Theses relationships represent a set of MxN homogeneous algebraic equations. In order 
to avoid a trivial solution (q^ = 0) , the nondimensional buckling load k is chosen such that 
the determinant of the coefficient matrix vanishes. The minimum value of k that satisfies the 
requirement is the buckling load. 


TABLE I. MATERIAL PROPERTIES 


Eu_ 


G„ 


9.25 Msi 

4.67 Msi 

2.27 Msi 

0.397 
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TABLE D. CONVERGENCE OF BUCKLING COEFFICIENT k FOR A SQUARE PLATE 


M/N 

1 

3 

5 

7 

1 

1.75 

1.69 

1.65 

1.65 

3 

1.72 

1.65 

1.65 

1.65 

5 

1.72 

1.65 

1.65 


7 

1.72 

1.65 



9 

1.72 





RESULTS 

First, a convergence study for Rayleigh-Ritz method is performed. Orthotropic plates with 
material properties as listed in Table I have been studied. For a square plate, buckling 
coefficients k versus increasing values of M and N are listed in Table II. It can be seen that 
when M=N=3, k reaches its final value. Convergence of the solution in this case is faster than in 
the case of long strips. For example a four-term series ( M = N = 3, note that only odd M and 
N are considered) converges to the final k value compared to a 12-term solution 
( M = 23, N = 1) for a long strip with /? = 20. 

The buckling coefficient k versus aspect ratio fi is plotted in Fig. 2. Figure 2 shows that 
for plates with aspect ratio less than six, the buckling coefficient k increases as /? increases. 
However, for long strips (/f > 6), k does not significantly change and takes a value of 
approximately S.75. 



Figure 2. Buckling coefficient versus aspect ratio 
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TABLE III. CONVERGENCE OF BUCKLING COEFFICIENT FOR A 
SQUARE PLATE USING GALERKIN METHOD. 


M 

k 

1 

2.54 

2 

1.96 

3 

1.79 

4 

1.68 

5 

1.66 


Galerkin method is also used to calculated the buckling load. The accuracy of the 
polynomial representation of the displacement field in the longitudinal direction, given by Eq 
(14), is examined for a square plate. The values of the buckling coefficient k are listed in Table 
UCconveiging to a value of 1.66 (M =5, N=l) compared to 1.65 as obtained using Fourier 
Series (Rayleigh-Ritz method, Table II). For M- 5 and N- 1, the displacement w, Eq. (14) 
takes the following form 

-V = [«, ,( 5 4 - 2« 3 + «) + 42 , ( i? 5 - 10$ 3 + 7 $)+ , 3 , - 12£ 5 + 5 { 3 ) 

+44l(4 7 -64*+ 5{ 3 ) + 45m( 4 5 8 -22{ 5 + 18^ 4 )]sin(OTj) (16) 


The buckling load intensity (load per unit width) calculated based on a parabolic and 
uniform load distribution is listed in Table IV. The parabolic load values represent the buckling 

load intensity at mid-span of the plate (i.e., ta„ or kl\ x it 2 l2b 2 ). 

>=a/2 

Table IV shows that the buckling load intensity is influenced by both the load distribution 
and the aspect ratio. For a square plate, the buckling load calculated based on a uniform stress 
distribution is 244% higher than the buckling load calculated based on a 
tensile/compressive/tensile parabolic load distribution. The difference is only 5% (reduction) 


TABLE IV. COMPARISON OF BUCKLING COEFFICIENTS FOR VARYING 
AXIAL STRESS AND UNIFORM STRESS 



Parabolic Load* 
(lb/in) 

Uniform 

Load* 

(lh/in) 

% diff. 


Fourier 

Series 

Polyno- 

mial 

Exact 




||g|| 

4.64 

+244% 

a=150 in. 
b=4.3 in. 
(a/b=35) 

m 

• •• 

5482 

-5% 


* Associated with shear and transvase stresses as given in Eqs. (4) and (5). 

* No shear and transverse stresses. 

* Buckling Load intensity at mid-span (x=a/2). ^ 
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for long strips. One may expect that a tensile stress at the edges as shown in Fig. 1 will tend to 
restore the plate to its stable position and therefore increases its buckling load. However, the 
parabolic stress distribution is associated with both compressive and shear stresses The 
transverses compressive stresses a y (y ) have a detrimental effect on the buckling load for 

square plates, while their effect is negligible for long strips. 


CONCLUSIONS 

The elastic buckling load of orthotropic plates subjected to a parabolically varying tension- 
compression-tension longitudinal stress, equilibrated by nonuniform shear stresses and 
transverse compressive stresses is analyzed. It is shown that the buckling load depends on the 
aspect ratio for plates with aspect ratio less than six. For square plates, the buckling load due to 
a uniform stress distribution is 244% higher than the buckling load calculated based on a 
“fensile/compressive/tensile parabolic load distribution. However, the difference is only 5% 
(reduction) for long strips. 
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